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Localized contractile configurations or asters spontaneously appear and disappear as emergent 
structures in the collective stochastic dynamics of active polar actomyosin filaments. Passive parti¬ 
cles which (un)bind to the active filaments get advected into the asters, forming transient clusters. 
We study the phase segregation of such passive advective scalars in a medium of dynamic asters, 
as a function of the aster density and the ratio of the rates of aster remodeling to particle dif¬ 
fusion. The dynamics of coarsening shows a violation of Porod behaviour, the growing domains 
have diffuse interfaces and low interfacial tension. The phase segregated steady state shows strong 
macroscopic fluctuations characterized by multiscaling and intermittency, signifying rapid reorga¬ 
nization of macroscopic structures. We expect these unique nonequilibrium features to manifest in 
the actin-dependent molecular clustering at the cell surface. 


Collections of active particles exhibit large concentra¬ 
tion fluctuations, leading in many cases to clumping and 
segregation into high and low density phases [1, 2] - the 
dynamics towards segregation and the space-time fluc¬ 
tuations of the order parameter, in single [3-7] and two 
component [8, 9] active systems, are found to be very 
different from their equilibrium counterparts. 

Recent studies on actomyosin-dependent molecular 
clustering at the cell surface [10-14] motivate us to ex¬ 
plore a novel kind of phase segregation driven by activ¬ 
ity. These studies suggest that the cortical layer of acto¬ 
myosin generates active stochastic stresses and currents 
at the cell membrane, driving local clustering of several 
membrane components [10, 12, 14]. The spatiotemporal 
statistics of such clustering, driven by an active noise, 
is well described using an active hydrodynamics frame¬ 
work [14], where membrane components that (un)bind to 
actomyosin are driven by the active filaments - we will 
call these components passive advective scalars or simply 
passive particles [14, 15]. The simultaneous presence of a 
dense stationary act in meshwork at the cortex [16] pro¬ 
vides a momentum sink to the dynamics of the cell mem¬ 
brane components. Membrane curvature effects appear 
unimportant, for such clustering occurs predominantly in 
the flat regions of the cell surface [10]. 

Here we study the potential scale-dependent segrega¬ 
tion of passive molecules from inert molecules (which do 
not bind to cortical actin) at the cell surface. Our model 
system consists of passive and inert particles embedded 
in an active medium in two dimensions (2d). While we 
observe both micro and macro phase segregation upon 
varying the active noise strength and active filament con¬ 
centration in the medium, the segregation dynamics and 
the nature of the active phase segregated state are dra¬ 
matically different from conventional phase segregation 
[17]. First, the active phase segregation is driven by ac¬ 
tive advection and not by gradients of chemical potential, 
and so occurs at temperatures higher than the equilib¬ 
rium segregation, a feature that is relevant to the cell 


surface context [18, 19]. Second, the domain coarsening 
dynamics shows strong departures from Porod behaviour 
[17], which implies that the interfacial widths are diffuse 
with low interfacial tension. Third, the active segregated 
phase shows strong fluctuations of the integrated order 
parameter and domain size, characterised by broad non- 
Gaussian distributions, in contrast to the sharply peaked 
distribution in closed systems. These strong fluctuations 
manifest as multiscaling and intermittency in the steady 
state, features seen in many driven nonequilibrium sys¬ 
tems such as turbulence [20]. As a result, the steady state 
exhibits a continual breakup and reformation of macro- 
scopically large structures. This fluctuation dominated 
phase ordering (FDPO ) is similar to those observed in a 
variety of semiautonomously coupled dynamical systems, 
such as passive sliders on a fluctuating surface [21], Burg¬ 
ers’ fluids [22] and active nematics [23]. 

Dynamics of active fluids - asters as emergent struc¬ 
tures - The basic hydrodynamic equations describing the 
polar contractile active fluid in terms of the active fila¬ 
ment concentration c, the polar orientation n, the apolar 
nematic orientation Q and the hydrodynamic velocity 
field v are described in recent reviews [2]. The precise 
form of the equations relevant to the present context are 
described in [14, 24]. Many studies, including [14, 24] 
have shown that the homogeneous, orientationally disor¬ 
dered configuration of active polar filaments is unstable 
to clumping and local orientation order to form spatially 
localized contractile regions - these emergent structures 
[25, 26] have been referred to as asters [14, 24]. 

Introducing noise parametrised by an “active tempera¬ 
ture” T a , results in a stochastic breakup and re-formation 
of asters (remodeling), with an exponential distribution 
of lifetimes [14, 15, 24] whose mean r r , is set by T a . This 
leads to a picture of a finite density of contractile struc¬ 
tures or asters, whose spatiotemporal fluctuations corre¬ 
spond to birth-death processes [15]. 

Passive particles that (un)bind to the filaments are 
driven by these contractile flows; when bound they are 
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FIG. 1. (a-b) Schematic of a binary system of passive A (blue) 
and inert B (chrome) particles in an active medium, (a) The 
active medium consists of polar active filaments (red arrows) 
and asters (circled), within the asters the arrows point in¬ 
wards, while outside the asters they are orientationally un¬ 
correlated. (b) The active medium consists of act in filaments 
(red) which form a mesh of actomyosin contractile platforms 
(circled), which draw in A particles that are bound to it. (c) 
PDF of the order parameter <f> for co — 0.128 and different 
values of A* at a fixed temperature T/T c = 1.17, where the 
equilibrium (Eqm) P{f>) shows a peak at <f> = 0. The bimodal 
distribution for oo > 0 indicates phase segregation, even at 
T > T c . (d) Phase diagram in (T, A*,xa) at oo — 0.128, 
shows micro (blue) and macro (red) phase segregation, when 
T > T c . 


advected with a current J oc cn; when unbound, they 
move diffusively. Thus the formation of contractile asters 
drive the localised clustering of passive particles (Fig. 
la,b). The stochastic remodeling of asters give rise to 
dynamical fluctuations of clusters of passive particles ex¬ 
plored in [14]. Here, we study the dynamics of phase 
segregation of passive-inert particles moving in an ac¬ 
tive medium consisting of a finite density of remodeling 
asters - however, rather than using the continuum hy¬ 
drodynamic equations [14, 24], we work with an effective 
theory that incorporates all its relevant features. 

Effective dynamics using kinetic Monte Carlo - Con¬ 
sider a collection of contractile platforms, represented as 
discs of diameter P*, distributed uniformly with area 
fraction A*, in a 2d system of size L. In one model of the 
active fluid medium, the filaments are distributed uni¬ 
formly and point radially inward within the aster, while 
they are oriented randomly outside the asters (Fig. la). 
Alternatively, our model could describe a situation where 
the active medium consists of patches of contractile ac¬ 
tomyosin mesh (Fig. lb). In either case, the contrac¬ 
tile platforms (e.g., asters) remodel with lifetimes chosen 
from an exponential distribution exp [—t/r r \. In this dy¬ 
namical background, consider a two component system 
of passive particles A (at concentration, x^) and inert 
particles B (concentration, x# = 1 — x^). The particles 
interact via a Lennard-Jones (LJ) short-range repulsion; 
in addition, A — A and B — B have a mutual attractive 
potential, u(r) = —J (a/r) 6 ,where J is in units of tem¬ 
perature ksT and a is the LJ-length scale. 

The passive particles that bind/unbind to the active 


filaments with rates k & (k u ), get advected with speed vq 
and timescale r v = R*/v o in a direction pointing radi¬ 
ally inwards, when the filaments are part of an aster. 
Outside the aster, the passive particle moves diffusively, 
with a diffusion coefficient set by vq and the orientational 
decorrelation time of the filaments. 

The kinetic Monte-Carlo (MC) dynamics [27] involve 
(i) Kawasaki exchanges [28, 29] of components A and 
B which obey detailed-balance, (ii) detailed-balance vio¬ 
lating active advective updates whenever A moves into 
an aster and (iii) detailed-balance violating stochastic 
birth-death updates of asters with a lifetime distribu¬ 
tion, exp [—t/r r \ [30]. One MC step is N- attempts at 
exchange, N being the total number of particles in the 
system. We have converted all time, length and energy 
units to s, /im and &bT, respectively [27]. We explore the 
dynamics as a function of particle concentration x^ and 
the active parameters characterizing the medium, namely 
oo = rf 1 /rf 1 , the ratio of the aster remodeling rate to 
the particle diffusion rate (jd is the diffusion time of par¬ 
ticles over the scale of an aster), aster area fraction A* 
and the Peclet number P e = rf 1 /rf 1 . We fix P e ^ 187 
and duty ratio K = k^/ + k u ) = 1 [27]. 

We find that the steady state distribution is very sen¬ 
sitive to both active parameters, oo and A*. At oo = 0, 
the steady state is equivalent to particles diffusing in a 
quenched random media of stationary asters, character¬ 
ized by an equilibrium distribution. At the other ex¬ 
treme, oo 1, the asters remodel much faster than the 
particle movement, the steady state is essentially that of 
equilibrium diffusing particles. It is only for 0 < oo 1, 
when the particles move fast compared to aster remodel¬ 
ing, that we get a genuine nonequilibrium steady state, 
where clusters form, fragment and move towards newly 
formed asters. 

Phase diagram - We first obtain the equilibrium phase 
diagram in the (T — x^) plane and the critical tempera¬ 
ture ksTjJ = 1.43, by monitoring the probability dis¬ 
tribution function (PDF) P(</>) of the coarse-grained seg¬ 
regation order parameter, (f> = , where Pa,Pb are 

the local densities of A and P, respectively. The order 
parameter cf> is coarse-grained over a scale 3cr which is 
smaller than the typical domain size [27, 28]. The phase 
boundary is determined by the appearance of a distinct 
double peak in P(</>) at (f> = ±1. 

We next determine the nonequilbrium phase diagram, 
as a function of (T, A*,x^) at different values of oo. Un¬ 
like the equilibrium case, phase segregation here is driven 
by the active advection of one of the components into 
asters, rather than by gradients in chemical potential, 
which implies that segregation can occur even at T > T c . 
Its extent depends on A* - for small A* only micro-phase 
segregation takes place (the domain sizes are roughly the 
aster size); as A* increases to a finite fraction of the sys¬ 
tem size, we start getting macro-phase segregation. Fig¬ 
ure lc,d indicates that the nature of segregation depends 
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FIG. 2. Active coarsening dynamics, shown here for uj — 
0.128, A* = 0.3 and T/T c — 1.17 : (a) Dynamical scaling of 
the order parameter correlation function g(r , t) with domain 
size R(t). The finite intercept of g(r, t) at the scaled r/R(t) —)> 
0 is an indication of long-range order, (b) R(t ) grows with 
time (red) and saturates to a value that depends on A* and uj 
[27], unlike the equilibrium phase segregation at T < T c (‘+’ 
symbols), where the domain size keeps growing as t 1 ^ 3 (black 
symbols in inset) until it reaches system size. In contrast, 
for active phase segregation, R(t) ~ t 1 ^ 4 before saturating 
(purple line), here uj = 0.128, A* as 0.6. (c) Structure factor 
S(k, t) also exhibits dynamical scaling but a violation of Porod 
behaviour - S(kR(t)) ~ /c _ ( 2 +o0 at large k, where 0 < a < 1. 
(d) a varies non-monotonically with uj, showing “equilibrium” 
behaviour (represented by the dashed line) at uj = 0 and 
uj » 1 . 


on the values of uj and A*. 

Coarsening dynamics - We next study the dynamics 
of segregation following a quench from the high temper¬ 
ature mixed phase to a region between the active and 
equilibrium phase boundaries (red and blue shaded re¬ 
gions in Fig. Id) - a situation relevant to the cell mem¬ 
brane, where equilibrium lipid phase segregation has only 
been observed at temperatures much lower than physi¬ 
ological temperatures [18, 19]. We study the dynam¬ 
ics of the order parameter correlation function g(r, t ) = 
(0(0, t)<j)(r , t)) as the system coarsens. As in conventional 
phase ordering [17], g(r,t ) appears to obey dynamical 
scaling (Fig. 2a, Fig. S2), g(r,t) = g(r/R(t )), where the 
domain size R(t) is obtained from the first-zero of g(r, t). 
The mean domain size R{t) grows before saturating to a 
value that depends on uj and A* - a power law fit yields 
(Fig. 2b, Fig. S3), R{t) ~ t 1 / 4 , different from the conven¬ 
tional phase ordering growth of t 1 / 3 [17]. Note that the 
maximal domain size is much smaller than the equilib¬ 
rium case, even though it increases with increasing sys¬ 
tem size L (with A* held fixed). The fourier transform 
of g(r,t), the structure factor S{k,t) also exhibits scal¬ 
ing S(k,t) = [R(t)] 2 S(kR(t)) (Fig. 2c, Fig. S4). In con¬ 
ventional phase segregation, S(kR(t )) at large k goes as 


k~ 3 , known as Porod’s law [17], signifying that the in¬ 
terface separating the growing phases is sharp. Here, in 
contrast, the tail of S(kR(t)) appears to go as k~( 2 + a \ 
where a depends on uj and lies in the range 0 < a < 1 
(Fig. 2d). This departure from Porod signifies that the 
interfaces are diffuse with very low inter facial tension, 
the interfacial width being £ ~ R 1_a (while still having a 
well-defined “domain”, i/R <C 1 for large R). This also 
implies that g(r/R(t)) has a cusp nonanalyticity at small 
r/R(t ), going as g(y) ~ 1 — y a ^ +_ 


Fluctuation dominated phase ordering - The most dra¬ 
matic aspect of the active segregated steady state is in the 
nature of fluctuations of a variety of statistical quantities. 
Supplemental Movies [27] show that the dynamics of 0- 
profile is characterized by strong fluctuations both at the 
interface and bulk, unlike the equilibrium situation where 
the fluctuations are appreciable only at the interface. To 
get a quantitative handle, we study the time series of R, 
the region-wise order parameter, 4/ = f Q 0(r) dr (where 
the integration is done over a large enough region of size 
Q <C L), and the fourier amplitude of 0 at low k (Fig. S5- 
S7), $(1) = m = 2n/L)\, $(2) = \4>{k = 4tt/L)|, 
<f>(3) = 10(fc = 6tt/L)\, — We find that the correspond¬ 
ing PDFs at steady state, show large variance and heavy 
non-Gaussian tails (Fig. 3, Fig. S5). Figure3d shows 

that the PDF P(S^) ~ exp —QA(uj) , ex¬ 

pressed in the large deviation function form [31], where 
— (\h). These observations are consistent with 
other driven non-equilibrium systems, such as passive 
sliders on fluctuating interfaces [21] and Burgers’ fluids 
[ 22 ]. 

The time series of various statistical quantities X (t) al¬ 
ternates between periods of quiescence and large changes 
over very short time scales (Fig. 4a), indicating a break¬ 
down of self-similarity due to intermittency. We char¬ 
acterize this behaviour in terms of the n-point structure 
functions S n = ((X(t) — X(0)) n ), where (...) is an aver¬ 
age over histories, as in studies of turbulent systems [20]. 
Here we compute such structure functions for R(t). 


These structure functions are expected to scale with 
time as S n ~ t nin (Fig. S8). In equilibrium phase segre¬ 
gation (away from critical points), the statistical quanti¬ 
ties are self-similar, implying y n = const - called gap scal¬ 
ing (Fig. S8). In the active phase segregation, y n is a non¬ 
linear function of n - a phenomenon known as multiscal¬ 
ing (Fig. S8). This breakdown of self-similarity is a conse¬ 
quence of intermittency (which is dominated by extreme 
events), as seen by the behaviour of “flatness” or kurto- 
sis, ft(£) = S^{t)/S^t) (Fig. 4b). Intermittence shows up 
as a divergence in n as t/r 0 (Fig. 4b), where r is a 
time scale which characterizes the lifetime of the largest 
structures in the system [34]. We argue that r scales 
as the remodeling time r r for large r r and is set by 
for small r r - this suggests the following nonlinear form, 
r = r r [uj 2 + cjq] , where ujq is a numerical fit parameter 
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FIG. 3. Time series of statistical quantities show large fluc¬ 
tuations in the phase segregated steady state, (a) Variance 
of R(t), A R(t) = ((R — (R)) 2 ), shows a non-monotonic de¬ 
pendence on uj (open circles). Line through data points is a 
piece-wise cubic spline fit. For comparison we show A R{t) for 
the equilibrium segregation state when T <T C (filled square). 
(b) Time series of region-wise order parameter T, evaluated 
over the square patch of size £4 (inset), for equilibrium and 
active segregation states. Color-bar indicates the value of the 
local order parameter cj). The mean value of T is smaller and 
the fluctuations are larger in the active state, (c) PDF P(£T), 
where SAt = T — (T), at different uj’s for fixed A*, shows 
a much broader width than that for equilibrium segregated 
state, (d) P(SAt) shows scaling with the region-size Q, the 
collapse curve is a gaussian. Data displayed for T/T c = 1.17. 
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FIG. 4. Strong fluctuations in steady state characterized by 
intermittency and multiscaling, (a) Time series of domain 
size R for different values of u j and A* shows sudden pre¬ 
cipitous drops and rises, suggesting intermittent behaviour. 

(b) Kurtosis n{t) as a function of t, for different values of uj. 

(c) Scaling collapse of n(t) with scaling variable t/r, where 
t — r r (uj 2 + ujq ) and the numerical parameter ujo — 0.1 
[32] at A* = 0.6, shows divergence as t/r -A 0. (d) Region 
(chrome) in uj — A* plane which exhibits fluctuation domi¬ 
nated phase ordering (FDPO). There is no macroscopic phase 
ordering in the blue regions. Data displayed for T/T c — 1.17. 


[32] (Fig. 4c). The onset of divergence appears earlier for 
faster remodeling r r . Note that n{t) shows “equilibrium¬ 
like” non-intermittent behaviour when uj = 0 and for 
uj > 1 (Fig. 4c). This is in agreement with our earlier 
assertion that these limiting cases correspond to steady 
states of an equilibrium model. 

Non-Porod behaviour, non-Gaussian distributions, 
multiscaling and intermittency are signature features of a 
driven nonequilibrium steady state called the fluctuation- 
dominated phase ordered (FDPO) state [21]. The most 
dramatic manifestation of the FDPO, however, appears 
in the continual remodeling of the largest scale struc¬ 
tures at steady state ( Movie S2). This visual feature is 
captured by the time series of the dominant fourier am¬ 
plitudes of the order parameter 4>(ra), for m = 1, 2, 3,.... 
The breakup of the larger structure into 2 or 3 domains 
is reflected in the anticorrelation of the cross-correlation 
function (<F(1, 0)4>(2, t )) and (<F(1, 0)4>(3, t)) (Fig. S6, S7) 
- the decrease in the amplitude <F(1) feeds into the ampli¬ 
tudes 4>(2) and <F(3) [35]. In Fig. 4d, we depict the region 
in the steady state phase diagram in uj — A*, where FDPO 
obtains. 

Discussion - In this paper we have studied the dy¬ 
namics of phase segregation of passive and inert parti¬ 
cles in the background of a fluctuating active medium in 
2d. We find that the active driving can induce segrega¬ 
tion at T > T c , the equilibrium transition temperature. 
The dynamics of coarsening, described by domain growth 
and sharpness of growing interfaces is very distinct from 
its equilibrium counterpart. More dramatic is the na¬ 
ture of the phase segregated steady state, which shows 
macroscopic fluctuations of the order parameter, mul¬ 
tiscaling and intermittency, characteristics of a FDPO 
state (Fig. 4d). 

We have now extended this study to temperatures 
T < T c , where the system can undergo an equilibrium 
phase segregation in the absence of active driving. Here 
again, we find that turning on activity results in a break¬ 
down of conventional segregation dynamics - R(t) grows 
and saturates to a value determined by A* and uj and the 
segregated phase shows FDPO. Activity thus destroys 
the very large domains obtained in equilibrium phase 
segregation and makes them more dynamic and intermit¬ 
tent. While our motivation and focus here has been on 
the segregation of passive particles driven by actomyosin, 
the general features of our study should be applicable 
to other hybrid passive/active systems such as particles 
driven by motor-microtubule complexes [36]. We look 
forward to testing many of the predictions reported here 
in the scale-dependent actomyosin-based segregation of 
molecules at the cell surface [37] and in in-vitro realiza¬ 
tions of supported membranes attached to a thin active 
layer of actomyosin [38]. 

We thank S. Saha, D. Foster and S. Mayor, and es¬ 
pecially K. Vijaykumar, A. Dhar, M. Barma and D. Das 
for insightful discussions. 
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